Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN M?eX style file vl.4) 



Quasi-radial modes of rotating stars in general relativity 



Shin'ichirou Yoshida^ ^ and Yoshiharu Eriguchi^ 

^SISSA, Via Beirut 2-4, 34014 Trieste, Italy 

^Department of Earth Science and Astronomy, Graduate School of Arts and Sciences, University of Tokyo, 
Komaba, Meguro-ku, Tokyo 153-8902, Japan 



o 
o 
o 

(N 
Oh. 

00 
OO 



> 

vn 
m 

00 

o 

ON 

Oh: 

6 



Accepted, Received 



ABSTRACT 

By using the Cowling approximation, quasi-radial modes of rotating general relativis- 
tic stars are computed along equilibrium sequences from non-rotating to maximally 
rotating models. The eigenfrequencies of these modes are decreasing functions of the 
rotational frequency. The eigenfrequency curve of each mode as a function of the rota- 
tional frequency has discontinuities, which arise from the avoided crossing with other 
curves of axisymmetric modes. 
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1 INTRODUCTION 

Radial oscillations of non-rotating relativistic stars have 
been studied for over thirty years. Methods for obtaining 
their spectra have been well established (Bardeen, Thorne 
& Meltzer 1966 ; see also Chapter 26 of Misner et al. 1973 
and the references therein) , and have been applied to several 
equations of state (see for example Meltzer & Thorne 1966) . 
These works were mainly motivated by consideration of stel- 
lar stability because g eneral relativistic effects tend to de sta- 
bilize stellar models ( fowler 1964 ; Chandrasekhar 1964). 

On the other hand, the effect of rotation on stellar os- 
cillations is less well understood in a general relativistic 
context. As in the non-axisymmetric mode case, the slow 
rotation approximation has been the only accessible way 
f or investigating the eigen mode behaviour of rotating stars 
( Hartlc fc Friedmanl975 ). Recently, numerical relativistic 
hydrodynamic codes have been developed by several au- 
thors and some numerical simulations of rapidl y rot ating 
stars have b een c arried out. Stergioulas et al. (200C) and 
Font et al. ( ]200C| ) have shown that initial small perturba- 
tions around an equilibrium star evolved to a superposition 
of normal mode oscillations (Note that their hydrodynamic 
simulation is done in the fixed backg round spacetime. On 
the other hand, Shibata et al. (2000) have solved the full 
system of Einstein equations to investigate the dynamical 
stability of rapidly rotating stars.). 

Although the excitation and evolution of these modes in 
realistic situations should be investigated by time dependent 
hydrodynamic simulations, it is also important to study the 
mode behaviour along rotational equilibrium sequences by 
using linear perturbation theory. 

So far we have studied a few sets of non-axisymmetric 
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eigenmodes of rotating stars in general relativity (f- 
modes by Yoshida & Eriguchi 1997,1999 ; p-modes by 
Yoshida 1999 (unpublished; presented at the 9th Yukawa 
International Seminar 'BLACK HOLES AND GRAVITA- 
TIONAL WAVES - New Eyes in the 21st Century-') ). 
These results have been obtained within the Cowling ap- 
proximation in which Euler perturbations of the metric co- 
efficients have been neglected (see McDermott et al. 1983 
and Finn 1988 for a definition ; see Lindblom and Splin- 
ter 1990 for the accuracy of the method when applied to 
non-radial modes of spherical stars). Apart from some low 
order modes, these results are in good agreement with those 
of the full perturbation theory including metric perturba- 
tions (For comparison of the eigenfrequencies for slowly ro- 
tating stars, see Yoshida & Kojima 1997 ; for comparison 
of the neutral points of the CFS instability, see Yoshida & 
Eriguchi 1999 which compare the results with the one ob- 
tained by full computation of Stergioulas & Friedman 1998 
and Morsink et al. 1999) 

It is therefore natural to expect that the Cowling ap- 
proximation could also be successfully applied to the quasi- 
radial modes which are the smooth extensions of the radial 
modes of spherical stars to rotating stars. In the present pa- 
per, we study quasi-radial modes by using the Cowling ap- 
proximation. Contrary to the expectations, our results in- 
dicate that computations with this approximation can not 
reproduce the relativistic instability of spherical stars. This 
is plausible because the instability is essentially caused by 
the loss of balance between gravity and the pressure gradi- 
ent, and in calculations of it even the small corrections of 
gravity cannot be neglected. Moreover, the phase cancella- 
tion of the perturbed gravitational potentials, which may be 
effective in the case of non-axisymmetric modes, cannot be 
expected to happen for radial modes. 

See the Appendix of the present paper for the compar- 
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ison of two methods in the case of radial modes of non- 
rotating stars. 

Although the validity of Cowling approximation for ro- 
tating stars is not fully assessed, we here expect that at least 
a qualitative picture of the eigenmode dependence on stellar 
rotation could be studied by this approximation. 



2 RESULTS 

The equation of state used here is the polytropic one, 



p = itp 



£ = p + Np 



(1) 



where p, e and p are the rest mass density, energy density 
and pressure of the stellar matter, respectively. Geometrized 
units, c = G — 1, are adop ted in this paper as well as 
Mq = 1, following Font et al. (EoOOl). The constant iV is the 



polytropic index. The adiabatic exponent of the perturbed 
matter is assumed to coincide with 1 -I- The factor k is 
another constant. 

Each equilibrium sequence is computed with k and 
fixed. 

In the present study polar-like coordinates are used and 
the metric components are written as: 

ds^ = —e^'^dt^ + e^°' {dr^ + r^dO^) + e^^r^ siv? 6{d(f> — ujdt)^. 

(2) 

The rotational axis is located at sin 9 = 0. 

The coordinates used in the actual numerical computa- 
tion are surface-fitted ones {r*,0*) which are defined as: 



r*=r/Rs{e), 



(3) 



where r = Rs (6) is the form of the stellar surface in equilib- 
rium. 

The numerical method u sed h ere is basically the same as 
that in Yoshida & Eriguchi ( 1997 ) where non-axisymmetric 
modes were investigated. A minor modification is needed 
to obtain the axisymmetric modes. In the case of non- 
axisymmetric modes, the Eulerian variable 5p/{e + p) is ex- 
plicitly set to zero at the centre of the star {Sp is the Eule- 
rian variation of the pressure). In the case of axisymmetric 
modes, however, this is not the case since the regularity of 
the solution requires d{Sp)/dr to be zero at the stellar cen- 
tre. Therefore we simply modify the finite-difference scheme 
at the innermost grid points in our numerical code. More- 
over to avoid the coordinate singularity on the rotation axis, 
points on the axis are excluded from the computational re- 
gion. 

Most of the results shown in this paper are computed 
with a resolution of 40 uniformly distributed gridpoints in 
the radial r* direction and 10 in the angular 6* one. The 
computational region is a quarter of the meridional section 
of stars, thus the range of the radial and the angular coor- 
dinates are < r* < 1, < < tt/2. 

In Figs. |l| and |^ the eigenfrequencies of the axisymmet- 
ric modes are plotted against the rotational frequency of the 
equilibrium model. 

The model parameters are tabulated in Table ^ 

The eigenfrequency and the rotational frequency are 
normalized by -^/pT/Ztt, where pc is the central rest mass 
density of the models. The sequences 'F', 'Hi' and 'H2' are 
the fundamental, first and second overtones of the quasi- 
radial modes, respectively. 



Table 1. Parameters of the stellar model. Here pc is the rest mass 
density at the stellar centre, which is fixed as constant along the 
sequence. M and AI/R are the gravitational mass and the mass- 
to-radius ratio, where R is the circumferential radius. 
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Figure 1. Sequences of axisymmetric models. Furidamental(F) , 
first (Hi) and second {H2) overtones are plotted as well as some 
other axisymmetric modes. The polytropic parameters here are 
N = 0.5, K = 6.024 X 10", Pc = 1.781 x 10"^. Frequencies are nor- 
malized by pc/in. The rotational frequency frot at the mass- 
shedding limit of this sequence is 0.597 when measured in this 
unit. 



Some of axisymmetric f-modes and overtones of p- 
modes are also plotted. ^ These are the continuation of the 
corresponding f- and p-modes with I = L and m = where 
I and m are the indices of the ordinary spherical harmonics 
yim{(^,'-p)- The label refers to a mode corresponding to 
the pn-mode with degree I = L and order m = in the 



t Our numerical code assumes reflection symmetry of the eigen- 
modes with respect to the equatorial plane of the equilibrium 
star. As a result, modes with odd integer I cannot be computed. 
Considering the order of the even I modes in the flgure, sequences 
of ^/and may be located somewhere between the the corre- 
sponding p-modes with I = 2,4. 
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Figure 2. The same as Fig.^ except that the parameters for the 
sequence are N = 1.5, k = 4.349, pc = 8.1 X 10~*. The rota- 
tional frequency f^ot at the mass-shedding hmit of this sequence 
is 0.2923 when measured in this unit. 



non-rotating limit. Similarly, the mode with the label ^ f is 
the f-mode with / = L, m = in the non-rotating limit. 

Generally the quasi-radial mode sequence encounters 
other sequences of f- or p-modes. It is seen that the so-called 
avoided crossing occurs on these sequences (see Figj^ . This 
seems to be the general relativistic extension of what has 
been found for the oscillat ions o f rotating Newton ian st ars. 
See for example Clement (1986) and Unno et al. (1989) for 
the Newtonian cases. 

Fig.^ clearly shows the presence of the avoided crossing 
whereby two eigenfrequency curves approach smoothly, and 
then depart from each other without crossing. At the point 
of closest approach, the characteristics of the modes on each 
sequence exchange. Thus the sequence of a mode with given 
characteristics has a discontinuity there. 

There are several discontinuities along a sequence whose 
number depends on the mode order as well as on the 
polytropic parameters of the equilibrium star. For the se- 
lected poly tropic parameters, discontinuities appear in the 
rapidly rotating models whose rotational frequencies are 
over ~ 80 % of the mass-shedding limit. 



Figure 3. A close-up of Figjl| around the point of closest ap- 
proach of the Hi and *pi sequences. 



3.1 Convergence and accuracy 

Since we use a finite number of grid points, our results neces- 
sarily contain some errors. To see how much the results differ 
when the number of grid points is changed, we compute the 
fundamental quasi-radial mode sequence by using different 
grid numbers {M,N) = (40, 10), (40, 20), (80, 10) where M 
is the number of grid points in r*-direction and N is that 
in 6l*-direction. By extrapolating three eigenfrequencies ob- 
tained from these grid numbers, 1^(40, lo] i ^^[40, 20] i ^^[80,10] i we 
estimate the converged value of the eigenfrequency vo in the 
limit of infinitesimal mesh size to be: 



fO — 2(l^[40,201 + l'[80,10l) — 3i'[40,10] ■ 



(4) 



3 SOME REMARKS ON THE ANALYSIS 



The ratio of V[m,n] to uo can be treated as a measure of 
convergence of the results when the grid number is changed 
(Fig.^ . As can be seen from this figure, eigenfrequencies ob- 
tained by using (Af, A'') — (40, 10) mesh numbers (our stan- 
dard resolution) agree with the converged values to within 
3 percent. 

It is noted that the relative error of iy[so,io] is smaller 
than that of 1^(40, 20] in the slowly rotating cases, however for 
rapidly rotating cases, the error of !^[8o,io] becomes larger. 
This is because for the rapidly rotating cases, the deforma- 
tion of the stars from the spherical configuration is so large 
that the lack of angular resolution becomes the main source 
of inaccuracy. 

To check the accuracy of our two dimensional (2D) code. 
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Table 2. Eigenfrequencies of the lowest order quasi-radial modes. 
The equilibrium sequence is constructed with the equation of 
state p = Kp^^^^^ with polytropic parameters TV = 0.5, k = 
6.024 X 10*, Pc = 1-781 X 10~^. Here Oc is the baryon mass den- 
sity at the stellar centre (see Table hi). In the present case, the 
model is at the mass-shedding limit when f^ot = 0.597. The ro- 
tational frequency frot and the eigenfrequencies are normalized 



Table 3. Eigenfrequencies of the lowest order quasi-radial modes. 
Polytropic parameters are N = 1.5, k = 4.349, pc = 8.1 X 10~*. 
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we have compared the results obtained here with those ob- 
tained with the one dimensional (ID) code described in the 
Appendix (Table In the ID code we employ the standard 
scheme to solve the eigenvalue problem of the linear ordinary 
differential equation within the Cowling approximation. 

As seen from this table, two results agree well to within 
several percent. 



3.2 Eigenfunctions and classification of modes 

Comparing with the eigenfunctions of non-axisymmetric 
modes, we notice that the eigenfunctions of quasi-radial 
modes change their shapes significantly along rotational 
equilibrium sequences. For example, the radial distributions 
of the Eulerian pressure perturbation and the radial compo- 
nent of the velocity perturbation change considerably near 
the equatorial plane of the star. The number of radial nodes 
of these functions increases as the stellar rotation rate in- 
creases. On the other hand, the overall radial dependence 
of the ^-component of the velocity perturbation changes lit- 
tle as the star spins up. Therefore we can use the shape 
of this function as a tracer of the selected mode along the 
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Figure 4. Convergence of the code with increased numbers of 
grid points. The ratio of the eigenfrequency of the fundamental 
(F) mode computed with given grid numbers to the extrapolated 
frequency uo, yjyQ, is plotted against the rotational frequency 
normalized by the Kepler limit frequency. The polytropic pa- 
rameters of the equilibrium state are AT = 0.5 , K = lO'^ and 
Pc = 8.1 X 10-*. 
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Table 4. Eigenfrequencies of the three lowest order radial modes 
of non-rotating stars computed with the 2D Cowling code are 
compared with the results obtained with the ID Cowling code. 
The polytropic parameters are the same as in Fig.|l|. In the 2D 
computation, the numbers of radial and angular grid points are 
40 and 10, respectively. 



Results by ID code 
2D code 



F Hi 
2.009 3.482 
1.996 3.449 



H2 

4.917 
4.865 




Figure 5. Eigenfunctions of the Hi mode for a slowly rotat- 
ing configuration. The polytropic parameters are the same as in 
Fig.p The rotational frequency frot = 0.3626. Each curve shows 
the dependence of each eigenfunctions for a fixed valued of 6* 
The abscissa is the radial coordinate distance normalized by the 
equatorial radius. The curves with the largest r value at their 
right-hand end correspond to those in the equatorial plane. The 
functions shown are q = 5p/{p + €) (upper left) and three velocity 
components Vr (upper right), Ve (lower left) and (lower right), 
where the Eulerian perturbations are employed. The eigenfunc- 
tions are normalized so that Vr{r* /R = 1; 9* = 7r/2) = 1. 




Figure 6. Eigenfunctions of the *pi-mode for the same stellar 
model as in Fig. ^. 
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equilibrium sequence. The nomenclature of mode sequences 
is based upon the behaviour of the mode of a non-rotating 
star to which the sequence is continued smoothly. 



Figure 7. Eigenfunctions of the Hi mode for a rapidly rotating 
configuration. The rotational frequency of the equilibrium star is 
0.5489. 



3.3 Realistic neutron star models 

In addition to the polytropic case presented here, we have 
tried to compute quasi-radial modes of realistic neutron 
stars by using some of the candidate zero temperature equa- 
tions of state. However it was rather difficult to obtain full 
sequences of these modes with our method. As the star be- 
gins to rotate, the convergence to the quasi-radial modes 
suddenly becomes much more difficult. This may partly be 
due to the fact that these modes are sensitive (as compared 
with the non-radial modes) to the surface condition of the 
equilibrium star. Unfortunately, in the case of more realistic 
EOS, the adiabatic exponent is subject to large oscillations, 
becoming negative in some parts. These large variations de- 
crease considerably the accuracy of our method which then 
becomes inadequate. 
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dr 



(e + p){m + Aixpr*) 
r{r — 2m) ' 



dm ^ 2 
—— = 4nr e, 
dr 



dv 
dr 



m + 4TTpr 
r{r — 2m) 



(A3) 



(44) 



where m(r) is defined from e^'*' = (1 — 2m/r) ^ . 

The equation of state is assumed to be polytropic, p = 
np^'^Tf , where k and N are constants, and the rest-mass 
density p is related to e hy e — p + Np. 



Figure 8. Eigenfunctions of the ''pi-mode for the same stellar 
model as in Fig. ^. 
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APPENDIX A: COMPARISON OF RADIAL 
MODES IN THE FULL THEORY AND IN THE 
COWLING APPROXIMATION 

To test the validity and accuracy of the Cowling approxima- 
tion for a spherical configuration, we here present a compar- 
ison between the results obtained by the full perturbation 
theory and by the Cowling approximation for low order ra- 
dial modes. 



Al Equilibrium model 

The space-time of the equilibrium star is characterized by 
the following metric: 

ds"" = -e'"'^''^dt^ + e'^'^'dr^ + r'(d6»' + sin'' Odip^). [Al) 

The metric coefficients, stellar pressure p and energy 
density e are obtained by integrating the standard set of 



A2 Equations for radial oscillations 

For radial oscillations of spherical stars, only the Lagrangian 
displacement function ^(r) is needed to describe the physical 
perturbation (see Chapter 26 of Misner et al. 1973) . 

The equation of motion of the displacement (r^e""^ = 
Q in the full perturbation theory is: 



C" + AfC + BfC = 0, 



where Af and Bf are defined as 



Af 



P_ 

P 



and 



Bf 



e+ p 



(z.T+4- 



+ X' + 



BTrpe 



(A5) 



(46) 



(47) 



The prime after a variable refers its derivative with respect 
to r. Here a is the frequency and the adiabatic exponent F 
is defined by: 



r = 



e -I- p Ap 
p Ae' 



(48) 



where A represents the Lagrangian perturbation. In general 
the adiabatic exponent need not coincide with 1 -f 1/A'^. 

In the Cowling approximation, the equation of motion 
is instead (r^e^^S^ = rj): 



r?" -I- 4c77' + Ben = 0, 
where 4c and B^ are defined as: 



4c 



A + v , 



and 



e + p 



I II , 2 2 

[ — V + a e 



(49) 



(410) 



(411) 



These equations can be solved by the matching method: 
i.e., we have to search for which makes the Wronskian 
of the solutions, obtained by integrations from the stellar 
centre and from the surface, vanish at some matching point 
inside the star. 
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A3 Boundary conditions 



The boundary condition for tiie equation of oscillations at 
the centre of the star is regularity of the variables. It requires 
^,ri ^ as r ^ 0. 

At the surface of the star, we impose the boundary and 
regularity conditions which reduce to 



+ 4u r + a e \Q ■ 



0, 



^ + P / II , 2 2A-2I/X 



{A12) 



{A13) 



A4 Results 

Keeping the polytropic index N fixed, we vary the com- 
pactness (mass-to-radius ratio) of the model to construct an 
equilibrium sequence. 

In Figs. Al and A2 we show the typical sequences of the 
three lowest order modes. For all of the sequences, the eigen- 
frequency obtained by the Cowling approximation is larger 
than that from the full theory: i.e., the Cowling approxima- 
tion overestimates the stability of the star. Before the turning 
point,^ the two curves are nearly parallel. As expected, the 
relative error of the results obtained by the Cowling approx- 
imation becomes smaller for higher overtones. 
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Figure Al. Squared eigenfrequency cr^ normalized by MR~^ 
is plotted as a function of the compactness M/R of the star. 
Here M and R are the gravitational mass and the radius of the 
spherical star in the Schwarzschild coordinate, respectively. Solid 
lines denote the eigenfrequencies of radial modes in the full theory, 
whereas dashed lines are those for the Cowling approximation. 
The polytropic index N is 3/2 and the adiabatic exponent V is 
5/3. 
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Figure A2. The same as Fig. Al except that the polytropic 
index TV = 1/2 and adiabatic exponent F = 2. Note that F is not 
necessarily equal to 1 + 1/Af. 



■t Note that with this parametrization, the turning point does not 
correspond to the maximum mass configuration. Thus the zeroes 
of the fundamental radial modes are not the turning points. 
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